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Abstract 



It is shown that the 3-body trigonometric G% integrable system is exactly- 
solvable. If the configuration space is parametrized by certain symmetric func- 
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1. Introduction 

Consider the Hamiltonian describing a three-body system with pairwise and three-body 
interactions 



G 2 2 

fc= 



A' 



1 



, a (a* - x z ) 2 ' fc<; (x fc + ^-2x m ) 2 ^ 



where 5,51 are the coupling constants. This model, now known as the rational G2 model, 
was originally proposed by J. Wolfes IJ who studied its bound states and proved its exact 
solvability. The scattering problem for the same model was explored in Q. It is worth 
mentioning that in the limit g\ — ► 0, the model ([O]) becomes the celebrated Calogero 
model §. 

Using the Hamiltonian reduction method (see also the review P]) Olshanetsky and 
Perelomov have shown that the model ( |1.1| ) is nothing but a particular (rational) case of 
a more general model related to the root system of the exceptional group G2 ■ They called 
this model the G2 model and proved its complete integrability. However, the question of 
whether even a trigonometric version of the general G2 model 

u = -V &2 1 g V 1 91 V - 

° 2 " 2 ^ dz 2 4 ^ sin 2 (±(* fc - x,)) 4 sin 2 (±(^ + * z - 2x m )) ' (L2) 

is exactly-solvable remained open for many years. Recently, Quesne @ showed its exact 
solvability for the particular case g = 0. In the present paper we demonstrate the exact 



solvability of (1.2) for the general case g,g\ 7^ 0. We use a notion of exact solvability 
based on the existence of a flag of functional spaces with an explicit basis preserved by the 
Hamiltonian. A constructive criterion for exact-solvability consists on checking whether 
the flag is related to finite-dimensional representation spaces of a Lie algebra of differential 
operators 0. If this criterion is fulfilled, then the corresponding algebra is called the hidden 
algebra of the system studied, and therefore the Hamiltonian can be written in terms of 
the generators of this algebra. 

it was shown that the eigenfunctions of the iV-body Calogero and Sutherland mod- 
form the flag of finite-dimensional representation spaces of the algebra gl(N), 
realized by first order differential operators. The corresponding Hamiltonians were rewrit- 
ten as quadratic polynomials in the generators of the Borel subalgebra of gl(N), and 
the coupling constants appear only in the coefficients of these polynomials. Recently, it 
was shown that this statement holds for all the ABCD Olshanetsky-Perelomov integrable 
systems (and also for their SUSY generalizations which turn out to be associated to the 




hidden algebra gl(N\N - 1), see gOf). 

Here we show that for the 3-body Calogero-Sutherland models, there exists one more 
hidden algebra that we shall call g^ (see Appendix B), besides the above-mentioned 
hidden algebra gl(3). This property is maintained for the general rational G2 model (1.1) 
and also for the degenerate {g = 0) trigonometric G2 model. However, for the general 
trigonometric G2 case only the hidden algebra g^ remains, while gl (3) is no longer a 
hidden algebra. 
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The paper is organized as follows. In Section 2 the algebraic and, the gl(3) and g^ 
Lie-algebraic forms of the 3-body Calogero and Sutherland models are given. Section 
3 is devoted to the demonstration of a phenomenon, which we call 'complementarity', 
as a consequence of which the 3-body Calogero and Sutherland models are shown to be 
equivalent to degenerate G2 models corresponding to Ql.l| ) an d ( |1.2| ) with g = 0. The 
special coordinates leading to the algebraic and the g^ Lie-algebraic forms of a general 
G2 model are used in Section 4, in other to demonstrate the exact solvability of this 
model. A realization of the algebra gl(3) in terms of first order differential operators, 
acting on the plane, is given in Appendix A. Appendix B is devoted to a description of 
the infinite-dimensional algebra of differential operators g^ , admitting finite-dimensional 
representations in terms of polynomials in two variables. Finally, in Appendix C the 
explicit formulas for the first several eigenfunctions of the general G2 model are presented. 



2. Algebraic and Lie algebraic forms of the 3-body Calogero and 

Sutherland models 

In this section we represent the Hamiltonians of the 3-body Calogero and Sutherland 
models in an algebraic form, by making use of two different sets of translationally invariant, 
permutationally-symmetric coordinates. The procedure leads to two different Lie-algebraic 
forms: (i) in terms of the generators of the gl(3) algebra || and (ii) in terms of some 
generators of the infinite-dimensional algebra g@' (see Appendix B). 



2.1. Calogero model. The Hamiltonian of the 3-body Calogero model is defined by 

3 r o2 1 3 1 



ftCal = ^ Yl 
i=l 



9 2 , 2 2 



i<3 J 



where g = v(y — 1) > —\ is the coupling constant and uj is the harmonic oscillator 
frequency. The ground state eigenfunction is given by 

^\x) = A y (x)e-^ , (2.2) 

where A(x) = ]X<7 I 

the Vandermonde determinant and X 2 = X^ x ?- ^ was 
shown by Calogero that the eigenfunctions for this model can be expressed as 

¥(a;) = y { c \x)P c (x) , (2.3) 

where P c (x) is a polynomial symmetric under permutations of any two Xj's. The operator 
having these polynomials as eigenfunctions can be obtained by performing on (^^) the 
gauge rotation 

fcCal = -2(^ C) (x))- 1 H C al^ C) (x) . (2.4) 

In order to study the internal dynamics of the system we introduce the center-of-mass 
coordinate Y = Yl^=i x j an d the translation-invariant Perelomov relative coordinates 

Vi = Xi-^Y , i = 1,2,3 , (2.5) 



which obey the constraint y\ + y 2 + 2/3 = 0. To incorporate the permutation symmetry 
and the translational invariance we consider the coordinates ||: 

r 2 = ~y\ - 2/| - 2/12/2 , (2.6) 
t~3 = -2/12/2(2/1 +2/2) • (2.7) 

In terms of these coordinates we obtain, after extracting the center-of-mass motion, the 
first algebraic form for the Hamiltonian (| 

2 r 2<% 3 T 3 - l 4wT 2 + 2(1 + 3^)J6V 2 - 6WT36V3 . 

(2. 



/ical = -2r 2 ^ 2T2 - 6r 3 ^ 2T3 + -rfc^ - [4u,r 2 + 2(1 + 3z,)]d T2 - 6 W t 3 <9 T3 



It is worth emphasizing that for a fixed coupling constant g there are 2 different solutions 
for u, hence the operator Q2.8] ) has two different sets of polynomial eigenf unctions. In fact, 
as shown in ||, the operator Q2.8] ) can be rewritten in a Lie-algebraic form in terms of the 
g/(3)-algebra generators (see (A.l) where n = and x,y replaced by r 2 ,T3, respectively). 
The corresponding expression is 

2 

heal = -2J 2j2 J2 - 6J 3 3 J^" + ^</ 2 ,3^ 2 ,3 ~ 4a;J 2 2 - 2(1 + 3^)7^ - 6wJ 3j3 . 

6 (2.9) 



The polynomial eigenfunctions of (p.q)- (|2.9| ) have the form given in (A. 2) (see also Fig.l in 
Appendix A). Further note that the operator (|2.8| ) is invariant under T3 — > — T3, resulting 
from the reflection symmetry (x — > —x) of the original Hamiltonian ( |2.lD . We incorporate 
this symmetry by introducing the new coordinates 

Ai = r 2 , A 2 = r 3 2 , (2.10) 

and then the operator ( |2.8| ) takes another algebraic form given by 

frcal = -2A 1 ^ lAl -12A a ^ lAa + |A?A 2 ^ aAa (2.11) 

-[4u\i + 2(1 + 3u)]d Xl - {12co\ 2 - ^\\)d X2 . 

(cf.(2.8)). Now /ical can be rewritten in terms of the generators (B.2) of the algebra gl 2 x 
R 3 C g^ (see Appendix B with n = and x, y are replaced by Ai, A 2 , respectively) only, 
rather than in terms of gl(3). The corresponding Lie-algebraic form for the Hamiltonian 
is then given by 

h Ca \ = -2L 2 L l - UL'^L 1 + ^L 7 L 3 - 2(1 + ?>v)L l (2.12) 

-AujL 2 - 12uL 3 - -L 7 . 

3 



The polynomial eigenfunctions of (2.11)-(2.12) have the form (B.l) (see also Fig. 2 in 
Appendix B) 



Both Lie-algebraic forms (2.9), ( 2.12j ) do not contain the positive-grading generators 
J21J3 an d L 4 , respectively. This ensures the preservation of the corresponding flags 
of polynomial spaces (A. 2) and (B.l), and the exact solvability of the Calogero model 
following the criterion formulated in the Introduction. 



2.2. Sutherland model. The Hamiltonian for the 3-body Sutherland model is defined 
by 



9 ^2 T iZj 



uth 



1 



2 ^ dx, 

k=i k 



k<l sm 2 (±(x k -xi)) 



(2.13) 



where g = v{y — 1) > —4 is the coupling constant. The ground state of this Hamiltonian 



is 



^ Suth) (x) = (A^ tri9 \x)) u , 



(2.14) 



where ^ tri9 \x) = Ili<j I s ^ n \ { x i ~ x j)\ is the trigonometric analog of the Vandermonde 
determinant. Sutherland showed that for this model any eigenfunction can be written as 



(Suth) 



x)PJe l 



(2.15) 



where P s (e lx ) is a polynomial, symmetric under permutations of any two x^s. These 
polynomials are the so called Jack-Sutherland polynomials. As for the Calogero case, the 
operator which has these polynomials as eigenfunctions can be obtained by performing on 
( 2.13 ) a gauge transformation 

h Snth = -2(^ Suth) (x))- 1 H S uth^ Suth) (x) . (2.16) 

To exhibit the internal structure of this system, we introduce the translation-invariant, 
permutation-symmetric, periodic coordinates [|| 

i] 2 = -^[cos(ayi) + cos(ay 2 ) + cos(a(yi + y 2 )) - 3] , (2.17) 



a 
2 



-^[sin(ayi) + sin(ay 2 ) - sin(a(yi + y 2 ))] 



CY' 



(2.18) 



In the limit a — > 0, these coordinates become (2.6)— (2.7). After extracting the center-of- 
mass motion, the Hamiltonian ( [2.16| ) in these coordinates takes the algebraic form 

~ 2 ~ 4 4a 2 
2 ,z 24; '°'~ mm v " ' 3 

>,2 



h 



Suth 



(6 + ^-m)V3d^ 2V3 



(2.19) 



+(|% 2 " ^vl)dl V3 - [2(1 + 3v) + 2(1/ + \)a 2 m ]d m - 2{v + \)a 2 m d m 



a 

3 T L-V- ■ -> ■ -V ■ 3/- -J-'/. -v- ■ 3, 

It is worth emphasizing that similarly to the Calogero case, there are 2 different solutions 



for v for a fixed coupling constant g. Consequently the operator ( 2. IS ) has two different 
sets of polynomial eigenfunctions. As shown in ||, the operator ( 2. IS ) can be rewritten 
in a Lie-algebraic form in terms of the y/(3)-algebra generators (see (A.l) with n = and 
x,y are replaced by 772,773, respectively). This procedure yields 



^Suth 



"2^2,2^2 — 6^3,3^2 + ~zJo.zJ? 



a 



2,3^2.3 



2(l + 3u)J- + ^-J,° 9 J.° 



24 



3,2^3,2 



(2.20) 



-Q 



2^2,2^2,2 + 3^2,2^3,3 + 2^3,3^3,3 + ^(^ + Y2~)(^2,2 + ^3,3) 



The polynomial eigenfunctions of ( 2,19| )-( |2.20| ) have the form (A. 2) (see also Fig.l in 
Appendix A). 
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Also, as in the Calogero case, the Sutherland Hamiltonian possesses a reflection sym- 
metry (x — > —x). We can incorporate this symmetry in a new set of variables by defining 

0"1 = 772 , cr 2 = r]l . (2.21) 

In these variables the operator ( |2,16| ) becomes 



h 



Suth 



2 4 

-(2a! + —a 2 - —a 2 )d ( 



,2 

CT1CT1 



8cr . „n 

(12 + — ai )a 2 d 2 a 



(2.22) 



/° 2 
+(3^1(^2 



2a 2 a 2 )& i 



0"20"2 



[2(l + 3i/) + 2(^ + i)a 2 ai]a a 



+ 



4 ^ 



7 

(- + 4z^)a 2 o- 2 



which is another algebraic form for the Sutherland Hamiltonian (cf.(2.19)). The operator 
(2.22) can be represented in terms of the generators of the algebra g&> (see (B.2), (B.4) 
with n = and x, y are replaced by a\, a 2 , respectively): 

,L 2 L 2 8L 3 L 2 



hs 



nth 



-21? i} 



+—T 

24 



12I?L l + -L'L A -a z { 



7r3 



+ 



2(l + 3i/)L 1 - -L 7 



1 



2(z, + -)a 2 L 2 



+ 2L 6 L 6 ) 
4(, + l)a 2 L3 



(2.23) 



which is the g^ Lie-algebraic form of the Sutherland Hamiltonian. The polynomial eigen- 
functions of (|2.22 ) or ( |2.23| ) have the form (B.l) (see also Fig. 2 in Appendix B). The 
operator ( |2.23| ) contains now the generator T, given by (B.4), which does not belong to 
gl 2 ix i? 3 . Note that this generator T is not of a positive-grading. 

Since both Lie-algebraic forms fl2.20| ), ( 2.23| ) do not contain the positive-grading gen- 
erators J2 , J% and L 4 , hence the corresponding flags of the polynomial spaces (A. 2) and 
(B.l) are preserved. This demonstrates the exact solvability of the Sutherland model fol- 
lowing the criterion formulated in the Introduction. Observe that if we put to = and 
a = 0, the algebraic and the Lie-algebraic forms of the Sutherland model reduce to the 
corresponding forms of the Calogero model. 

3. Complementarity 

Let us consider a quantum-mechanical 3-body system. There are two systems of 
translation-invariant relative coordinates y and y related by: 



yi 

m 
m 



yi 
y-2 
2/3 



y-2 

2/3 

in 



Let the Hamiltonian for this 3-body system be given by 

d 2 /a 2 d 2 d 2 



n 



-A 



dY 2 



D 



where 



Vi(y) 



\dyj dy\ 

3 



Vi{y) + V 2 (y) 



(3.1) 



(3.2) 



Vtivi) , V 2 (y) = ^2v 2 (yi) . 
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Here Y is the center-of-mass coordinate and A, B are some constants. Suppose the co- 
ordinates y in (3.2) are the Perelomov coordinates ( |2.5| ), then one can easily see that 



(3.3) 



A = 3, B = | and the coordinates y are nothing but the Jacobi coordinates 



Vi 



On the other hand, if we take the y's as the Jacobi coordinates, then in (3^2) A = 3, B = 2 
and the y's become the Perelomov coordinates. So beginning from some Hamiltonian 
with a definite choice of relative coordinates, we arrive, after the change of variables 
( |3.1| ), at a Hamiltonian where the potentials V\ and V2 are interchanged. Of course, the 
spectra remain the same under such an operation, which makes it possible to connect 
the Hamiltonians with different coupling constants for two- and three-body interactions. 
Thus, we can put in correspondence two distinct physical problems by regarding them as 
related by some kind of 'complementarity' (equivalence). It is important to note that this 
peculiar feature appears only for 3-body problems. 

In order to illustrate the above discussed complementarity, let us consider the Hamil- 
tonian (2.22). Introduce the new coordinates (by replacing y — > y in ( |2.21| )) 

cos(a(yi - y 2 )) + cos(a(y 2 - 2/3)) + cos(a(y 3 - y x )) - 



0"! 



02 



0-2 {f) 



Q 

4 

a 6 



sm(a(yi - y 2 )) + sin(a(y 2 - 2/3)) + sin(a(j/ 3 - yi)) 



It can be verified that these coordinates are related to the original a coordinates 
by the following non-trivial algebraic relations 



(3.4) 
(3.5) 



<7i = 3(Ji + \a 2 a\ + \a A a 2 , 
2 o 

&2 = — 27(72 — 4(7]* — 9a 2 0"i<72 



24 I42 1 6 2 

a a 1 - -a a x a 2 - z^a D a 2 ■ 
2 lb 



(3.6) 
(3.7) 



Adding back the center-of-mass coordinate to the Hamiltonian (2.22) with ui, CT2 replaced 
by <5"i , 02 , followed by the gauge transformation 



H 3B = ^ B) (x)h Suth (a 1 ,a2)(^ D '^)) 



(3B), 



where the gauge factor is defined by 



(3B) 




(x) = (Af ri9 \x)) 



(3.8) 



(3.9) 



with A\ u iy> (x) = n?<j i j^k I srn \ ( x i + x j ~ 2xk)\ , we finally arrive at the 3-body Hamil- 



tonian with the 3-body interaction studied by Quesne 

2 



^3B 



k=l 



dxl 4 



3 

E 

k<l 
k,lj^m 



sin (oixk + xi- 2x m )) 



(3.10) 



Note that now the coupling constant g\ is given as 

9i = 3v(v - 1) . 

Thus, using the 'complementarity', we have been able to relate the 3-body Sutherland 
model with two-body interactions and the coupling constant g = v{y — 1) to a three-body 
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model with 3-body interactions and the coupling constant g\ = ?>v{ii — 1). Although the 
solvability of the last problem was already proved in || by finding the hidden algebra 
gl(3), the complementarity gives an immediate explanation for the solvability and a very 
short and transparent way of proving it. It also clarifies why the polynomial eigenfunctions 
of this model are given by Jack-Sutherland polynomials as noted in ||. 

4. Solvability of the G 2 integrable system 

Making use of the previously introduced relative coordinates, we shall derive in this 
section the algebraic and Lie-algebraic forms of the rational and trigonometric G2 models 
which lead to polynomial eigenfunctions. This fact together with the explicit calculation 
of the eigenvalues, exhibit the exact solvability of the model. 

4.1. The rational G2 model. As pointed out in the Introduction, this model was solved 
exactly for the bound states by J. Wolfes Q, who used the procedure of separation of 
variables. In this section the algebraic and Lie-algebraic nature of the solvability of this 
model will be shown. 

We begin by recalling the Hamiltonian for the rational G2 model 



k=l 



3 3 ^ 



( r ) 1 \ ^ O 09 \ 1 \ ^ 



~k<i (^-^) 2 f^l (Xk + x l -2x m )* ' (41) 

k,lj^m 

where g = v(v — 1) > — | and g\ = 3/i(/x — 1) > — | are the coupling constants associated 
with the 2-body and 3-body interactions, respectively. In this case we have two coupling 
constants and two different solutions for each value of v and /i, which implies an existence 
of four families of solutions. The ground-state eigenfunction is given by 

^ r) (x) = (A^(x)y(A { ;\x)Ye-^^^ , (4.2) 

where A^ r \x) = Ili<j \ x, i ~ x j\ an d A^[\x) = W^j. ij^k l x * + x j ~ ^ x k\- Wolfes showed 
explicitly that any solution to ([O]) can be written in factorizable form as 

*(x) = ^ (*)pg(x) , (4.3) 

where Pq^(x) is a polynomial symmetric under permutations of any two of the Xj's. The 
operator having these polynomials as eigenfunctions can be obtained by gauge rotating 
(0): 

/ l g = -2(^(x))- 1 « ) (.). (4.4) 

It then follows that in terms of the coordinates r, given in (|2]D-(|3), the operator flP|) 
takes the form 

h%\ = -2r 2 d 2 T2T2 - 6r 3 d 2 T2T3 + |rf <% T3 - {4wr 2 + 2[1 + 3(fi + v)]}d T2 - 6^r 3 5 r3 , 

3 (4.5) 



after dropping out the center-of-mass dependence (cf. (2.8)). This operator is the first 



algebraic form of the rational G2 model. It can be immediately rewritten in terms of the 
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generators of the algebra gl(3) (see (A.l) where n = and x,y are replaced by t 2 ,T3, 
respectively) 

2 

!'■' !'■■> i . /'■' 'H I i *.>/,. i . , \ I / /,' . I ' 

1,3 ? 

(4.6) 



h ( G ] 2 = -2J 2 ° 2 J 2 " - 6J 3 ° )3 J 2 - + - J 2 ° 3 J 2 ° 3 " W 2 ° 2 - 2[1 + 3(/x + !/)] J 2 - - 6a;J 3 °< 



(cf. ( p.9| ) ) . Eq. ( |4.6| ) is the 5/(3) Lie-algebraic form of the rational G 2 model. 

The model Kq admits another set of algebraic and Lie-algebraic forms. They can 
easily be obtained if the operator ( f4.4p is written in terms of the coordinates Ai, A 2 given 
by pT0|) . Thus, 

= -2A 1 ^ lAl -12A 2 ^ lA2 + ~A?A 2 df 2A2 (4.7) 



-|4wAi + 2[l + 30* + i/)]|9 Al - (!2wA 2 -^A?)9 A2 , 



(cf. fl2.11| )). This is the second algebraic form of the rational G 2 model, which also admits 
a representation in terms of the generators of the algebra g^ given in (B.2), where n = 
and x,y are replaced by Ai,A 2 , respectively, as 



h$ = -2L 2 L 1 - 12L 3 L X + -L 7 L 3 - 2[1 + 3(/i + v^L 1 - 4uL 2 - 12uL 3 - -L 7 

3 3 (4.8) 



(cf. ( p. 12 )). Eq. ( [4.8| ) is the g^ Lie-algebraic form of the rational G 2 model. The operator 
(4.8) depends on the generators of the gl 2 x i? 3 -algebra only. This implies that (4.8) 
possesses two invariant subspaces, W n and W n (see (B.l), (B.3)). Thus, it leads to the 
conclusion that there exists a family of eigenfunctions depending only on the variable Ai. 
This, in fact, is already known both for this model Q as well as for the general many-body 
Calogero model ||. This fact was used in [12| to construct quasi-exactly-solvable many- 



body problems which generalize the Calogero model. We should emphasize that the same 
property allows one to construct also quasi-exactly-solvable generalizations of the rational 
G 2 model. Such generalizations will be presented elsewhere. 

In conclusion, one can state that the rational G 2 model, as it is for the case of the 
Calogero and Sutherland models, admits two different algebraic and also two different Lie- 
algebraic forms and, eventually, is characterized by two different hidden algebras. Since 
the Lie-algebraic forms ( f4.6[) , (|4.8j ) contain no positive-grading generators, four infinite 
families of polynomial eigenfunctions of (fO|) and (|4.7| ) occur according to the different 
values of z/, fi. 

4.2. The trigonometric G 2 model. The Hamiltonian for the trigonometric G 2 model 
has the form 

- d 2 ,g 1 , gi 1 



n, 1 \ - a g \ - 1 gi \ ~> 



2 t^i dx l ' 4 ^sin 2 (i(x fc -^)) ' 4 ^ s\n 2 {\{x k + xi-2x m )) ' ^ 

where g = v(y — 1) > —\ and g\ = 3yu(/i — 1) > — | are the coupling constants associated 
with the 2-body and 3-body interactions, respectively. Again, as in the previous case, 
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there are four families of solutions depending on the permitted values for v and /i. The 
ground-state eigenfunction is given by 

V$\x) = {^ tri9 \x)f \&f' i9] \x)Y , (4.10) 

where A (:c), A? ri9) (x) are the trigonometric analogies of the Vandermonde determi- 
nant and are defined by 

A^\x) = Y[\ S m\(x i -x 3 )\ , 

i<j 



3 

Af rt9 \x)= Yl \ S m-(x i + x J -2x k )\ . 



k<l 

Guided by what occurs with the Calogero, Sutherland, rational and degenerate trigono- 
metric G 2 models, let us check whether there exist families of eigenfunctions of (|4.9|) of the 
type ^f(x) = VPq • > (x)P(g 2 )(x), where the P(g 2 ) are polynomials. If this is the case, there is 
a chance that following the conjectures made in Q the trigonometric G 2 model also pos- 
sesses a hidden algebraic structure. For this purpose we make the gauge transformation 
of ( [4.9D with the ground state eigenfunction Q4.1C| ) as the gauge factor. We get 

h Ga = -l^S\x))- 1 7U h (9§\x)). (4.11) 

Rewriting ( [4.11 ) in terms of the new coordinates a given in (|3.4|) -([3l]|) we arrive, after 
some calculations accompanied by some miraculous cancellations, at an operator with 
coefficients which are surprisingly polynomials. The resulting expression is 

h G2 = (2a x + — or?-— or 2 )9?. +(12 + — a^dj^ (4.12) 



-(~ala 2 - 2a 2 d 2 2 )dl 2 - a2 + {2[1 + 3(p + 2u)] + ~(1 + 3^ + Av^^jdg 

4 7 1 

-(1 + Au)a\ -[- + 4(/i + u)]a 2 a 2 \d& 2 , 



3 V ' 1 L 3 

(cf. fl2.11|) , fl2.22|) , (|4.7|)). Eq. ( 4.1 2| ) is the algebraic form of the trigonometric G 2 model, 
and can be written in terms of the generators of the algebra containing the generators 
of its subalgebra gl 2 x R 3 , plus the extra generator T. The explicit expression is given by 

h G2 = 2L 2 L 1 + 12L 3 L 1 --L 7 L 3 + a 2 (— + ^^ + 2L 3 L 3 )-^-T (4.13) 

3 2 3 24 

+2[1 + 3(a* + 2u)]L 1 + Au)L 7 

+ (2pi+~v+±)a 2 L 2 + [^ + 4(fi + ^]a 2 L 3 , 
where in the generators (B.2),(B.4) the parameter n = and x, y are replaced by efi, a 2 , re- 



spectively (cf.( |2.23| ),( |4.S| )). This is the Lie-algebraic form of the trigonometric G 2 model. 
However, unlike all previously discussed examples including the rational G 2 model, there 
does not exist anymore a gl(3) Lie-algebraic representation for the trigonometric G 2 model. 
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A partial explanation for this fact is related to the non-existence of the algebraic repre- 
sentation of £71,2 in terms of 012, similar to the representation (|3.6| )-( ^7| ). 

The Lie-algebraic form ( 4.13j ) contains no positive-grading generators such as L 4 in the 
algebra and, hence, the flag of polynomial spaces (B.l) is preserved. This demonstrates 
the exact solvability of the trigonometric G 2 model, once again following the criterion 
formulated in the Introduction. Note that there are four infinite families of polynomial 
eigenfunctions labeled by the different values of p. These polynomials are generalizations 
of the Jack-Sutherland polynomials. 

The next question is to find the spectrum of the Hamiltonian ( |4.9| ). First, let us note 
that the operator ( [4.12 ) has a block triangular form in the a variables. There is a simple 
trick allowing us to reduce this operator to pure triangular form, based on the introduction 
of the new variables 

Pi = d-i , P2 = o 2 + \o\ . (4.14) 

cr 

It is worth noticing that if a — > this change of variables becomes singular, reflecting 
the non-existence of bound states for the rational Calogero and G 2 models in the absence 
of the harmonic oscillator term in potential. This coordinate transformation has the very 
attractive property of leaving the space (B.l) invariant. 

In these new coordinates the Hamiltonian ( [4,12| ) takes the form 

h G2 = (2p 1 + ^a 2 p 2 -^p 2 )d 2 pipi + (l2p 2 + 2a 2 p 1 p 2 -^p 2 )d 2 pip2 (4.15) 

96 256 2 

+{2a 2 p 2 2 + —p lP2 - — p \)d 2 p2p2 + [2(1 + 3^ + 61/) + -(1 + 3/x + 4v)a 2 Pl ]d pi 

16 

+{2(1 + 2p + 2v)a 2 p 2 + -^(2 + 3p + 6v)pi}d P2 , 

or 

and it is easy to check that it is indeed a triangular operator. Evidently, this operator can 
be rewritten in terms of the g( 2 ' -generators as 

h G2 = 2L 2 L 1 + Y2L 1 !} - + ^a 2 L 2 L 2 + 2a 2 L 2 L 3 + 2a 2 L 3 L 3 (4.16) 

l| L 7 L i + ^ L 3 L 6 _ ^ L 6 L 7 + 2(1 + 3 + 6v)L i 
or or or 
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+(2 P + -v)o?L 2 + 4(/x + v)a 2 L 3 + ^(2 + 3/i + 6u)L b 
3 a A 

where in the generators (B.2),(B.4) the parameter n = and x,y are replaced by pi,p 2 , 
respectively. 

Using either one of the representations ( [4.15 ) or ( [4.16| ) the energy levels of the Hamil- 
tonian Hq 2 can be easily found and are given by 



where n, m are quantum numbers with 



2 4 2 

— (n — 2m — l)(n + m + 1) + 2np + ~^{2n — m)v + — 



o? 



(4.17) 



n = 0, 1,2,3,... , 0<m< Q) 
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The explicit expressions for the first several eigenfunctions of ( 4.1 5| ) are presented in Ap- 
pendix C. 

5. Conclusion 

We have found that the general trigonometric G2 integrable model with two arbitrary 
coupling constants is exactly-solvable. This model is characterized by a certain hidden 
algebra, which is an infinite-dimensional Lie algebra of the vector fields. It is quite amazing 
that for all degenerations of the general G2 model and also generalizations of the rational 
cases obtained by adding the harmonic oscillator interaction, there exists an alternative 
hidden algebra gl(3), which disappears only in the general case of two non- vanishing 
coupling constants. For completeness of the presentation, the explicit forms for some of 
the eigenvalues of the general trigonometric G2 integrable model are given. 

It is important to emphasize that the algebraic forms for the Hamiltonians of the 
Calogero, Sutherland, rational and trigonometric G2 integrable models (2.8), (2.11), (2.19), 
(2.22), (4.5), (4.7), (4.12), (4.15) are not contained in the list of algebraic forms of 
the Schroedinger operators possessing hidden algebraic structures, given in the papers 

mm mi. 
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Appendix A. Representation of the algebra gl(3) 

The algebra gl(3) has a realization in terms of first order differential operators acting 
on the (x, y)-plane which are given by 

J 2 - = a x (-1,0), J%, 2 = xd x (0,0), J$,2 = yd x (-i,+i), 

J 3 " = 8y (0, -1) , J 3 3 = yOy (0, 0) , J 2 3 = X8y ( + 1, ~1) , 

J° = n -xd x - yd y (0, 0) , J+ = xJ° (+1, 0) , J+ = y J° (+1, 0) , 

(A.l) 

where n € R. In this realization a grading (a, (3) can assigned to the generators (A.l) by 
observing that Jx p y q oc x p+a y q+/3 . If n is a non-negative integer, the representation (A.l) 
becomes finite-dimensional and the representation space 

V n = (xPyi\0<(p + q)<n), (A.2) 

is characterized by the Newton diagram (see Fig.l). 




n 



Fig.l The Newton diagram illustrating the finite-dimensional representation space (A.2) 

of the algebra gl (3). 

The algebra (A.l) acts on the space (A.2) irreducibly. 

Appendix B. Representation of the algebra g^> 

We define the algebra g^ as an infinite-dimensional algebra of differential operators on 
the (x, y)-plane possessing a finite-dimensional representation which acts irreducibly on 
the space of inhomogeneous polynomials 

W n = {xPyi\0<(p + 2q)<n}, (B.l) 

which is diagrammatically represented by Fig. 2. 
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Fig. 2 The Newton diagram illustrating the finite-dimensional representation space (B.l) 

of the algebras g^ and gl<i x B? . 
The algebra contains the subalgebra gl2 x B? of first order differential operators 



L 



Or. 



L 2 



cd T , 



n 
3 



(0,0) 



L 



L< 



yd y -- 

8y 

x 2 d v 



x 2 d x + 2xyd y 



nx 



(+1,0) 



xd„ 



(+1,-1) 



(-1,0) , 
(0,0) 

(0,-1) , 

-v (+2,-1), (B.2) 

where (,) denotes the grading (defined in the similar way as in Appendix A) and n is a 
non-negative integer. The gh X -R 3 -algebra, (B.2) is characterized by the commutation 
relations 

[L\I?\=L\ 
[L 2 ,L 3 ] = 0, 
[L 5+k , L 1 ] = -kL 5+k ^ , 
[L 5+k ,L 4 ] = (2-k)L 5+k 



[L\L*] = 

\l 5+ Kl 2 } 



o , 



kL 5 ^ 



[L 1 ,L 4 ] = 2(L 2 + L 3 ) , 
[L 3 ,L 4 ] = 0, 
[L 5+fe ,L 3 ] =0 , 
[L 5+k ,L 5+m ] = , 

and the generator L 4 is the only positive-root generator of the algebra gl2 x B? . If this 
generator is omitted the remaining generators form the Borel subalgebra of gl2 x i? 3 . The 
algebra gl2 x B? acts on (B.l) reducibly possessing an invariant subspace 

W n = (x p |0 <p < n) . (B.3) 

In order to obtain the Lie-algebraic forms of Calogero and rational GVintegrable models 
it is sufficient to use the gl2 x B? generators only representing their Hamiltonians. However, 
the Lie-algebraic forms of the Sutherland and trigonometric GVintegrable models require 
more generators than these the gl2 x B? algebra can provide. One of such extra generators 
leaving W n invariant has the form 

T = yd 2 xx , (B.4) 

and is characterized by the grading (— 2,+l). Clearly it does not belong to the universal 
enveloping algebra of the gl2 x i? 3 -algebra. The generator (B.4) does act on the space 
(|B.3|). Thus, the space (B.3) is no longer a common invariant subspace with respect to 



4 We are grateful to R. Bautista for fruitful discussion on the 3/2 X i? 3 -algebra 
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the action of the operators (B.2), flB.4| ), Consequently, the operators (B.2), (BA) act on 
(B.l) irreducibly and one can state that the g^ -algebra coincides with the algebra of all 
polynomials in the generators (B.2), \B.3Q . The proof is based on Burnside theorem (see 
discussion in |?]]). 

Due to the fact that the Hamiltonians we study are represented by second order differ- 
ential operators, it would be instructive to classify all second order differential operators 
acting on W n , or, equivalently, these having W n as the invariant subspace. It is clear 
that some of these operators are given by quadratic elements of the universal enveloping 
algebra of gl<i x -R 3 being quadratic polynomials in the generators (B.2). In order to find 
the remaining ones we should consider all possible the commutators, double-commutators 
etc between the operator T and the generators of the algebra gl% x R 3 . 

Firstly, let us check the commutation relations of T with the gfo x .Regenerators: 



\L\T] 



[L\T] 



, 

-2xyd: 



,2 

XX 



[L 5+k ,T] 



[L 2 ,T] = -2T , 
Ay 2 d 2 xy + 2(n - l)yd x 

k(k - l)x k - 2 yd y 



L 2 L l 
L 2 L 2 



2L 1 L 3 
4L 2 L 3 



[L 3 ,T] =T , 
T 1 (-1,+1) 

2kx k - l yd 2 xy 

k = 
k = 1 



(B.5) 



L 2 



(2 + |n)L 3 -(| + f)n , k = 2 . 



We find that except for the commutator [L 4 , T] =T l all other commutators belong to the 
universal enveloping algebra of 5/2 x i? 3 . Now among the double commutators we should 
consider only those involving the commutators of T 1 with the 5^2 x i? 3 -generators, namely 



[L\T 1 ] = 


-2T , 


[L\T 1 ] = 


-2L\L 2 + 4L 3 ) , 


[L 2 ,T 1 } = 


-T\ 


[L«,T 1 ] = 


-2L A L l + AL 3 L 3 - 2L 2 


[L\T 1 } = 


T\ 




2(1 + ^)L 3 - n(l + 2 -n) 


[T,T 1 } = 


, 


[L 7 ,T 1 ] = 


-2{L 2 - 2L 3 )L 4 + 2L 4 


[L\T l ] 


= 2x 2 yC 


+ 8xy 2 dl y 


+ 8 V 3d yy 



+ 4(n - l)xyd x + (12 - 8n)y 2 d y + 2n(n - l)y 



= T 2 (0,+l). 

Hence, the only new operator that appears is T 2 
we get 



(B.6) 



Repeating the above procedure with T 2 



[L\T 2 ] 
[L 3 ,T 2 ] 

[L 5 ,T 2 ] 
[L 6 ,T 2 ] 



-2T 1 



-2L l L A + 24L 3 L 3 + 12L 2 L 3 



2L 2 L 4 + 8L 3 L 4 



[L 2 ,T 2 ] 
[L\T 2 ] 

4(1 +n)L 3 
[L\T 2 \ 



, 

, 
2 

— n , 

3 ' 

2L 4 L 4 



[T , T 2 } 

[r 1 , r 2 ] 



, 
, 

(B.7) 
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As seen there are no new operators found. Therefore the three generators T, T 1 , T 2 exhaust 
all generators which can be represented by second order differential operators and which do 
not belong to the universal enveloping algebra of the gl2 x -R 3 algebra. We should emphasize 
the surprising fact that these new generators T, T 1 ,T 2 form an abelian subalgebra of 
(see (B.6), (B.7)). 
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Appendix C. Eigenfunctions of the G2 trigonometric model 

We present here the explicit expressions for the first six eigenfunctions and eigenvalues 
of the general trigonometric G2 model for the Hamiltonian ( [4.15 ). 



F ( 2 J_0 _L 8 ^1 2 

£i >0 = \-+2^+-u\a , 

1 + 3m + 61/ 

tin = pi H 

1 (1 + 3fi + 4i/)a 2 

( 8 16 \ j 

£2,0 = \ - + 4/J + —v I a 



2 a 2 3(4i/ + 1)(2 + 3m + 6i/) 9(4i/ + 1)(2 + 3/1 + 6i/)(l + 3/1 + &v) 

Pl 8 + \6v P2 + (1 + 2~vj(3 + 3/i + 4i/)a 2 Pl + (1 + 2i/)(3 + 3m + 4i/)(2 + 3m + 4i/)a 4 



£2,1 = ( - 2 + 4p + 4i/ I a 



2 a 2 13(2 + 3m + 61/) 13(2 + 3m + 6i/)(l + 3m + 61/) 

* 2,1 Pl + 48 P2 + 2(-4 + 3m + 2i^) P1 + 2(-l + 2m + 2i/)(-4 + 2v + 3/x)a 4 



£3,0 = ^6 + 6M + 8^a 2 
, 3a 2 

16(1 + I/) P II / I I ■ 1 -r : ■>/'.'<'- MM I. /t y —''.'!■■' >/' .1 ' ' 

9(1 + 4i/)(5 + 3m + 6i/)(2 + 3m + 6i/)(4 + 3m + 61/) 9(1 + 4i/)(5 + 3m + 6i/)(2 + 3m + 6i/)(4 + 3m + 6i/)(l + 3m + 61/) 

(1 + i/)(4 + 3m + 4!/)(2 + 2i/ + m)(5 + 4i/ + 3m)c* 4 Pl (1 + i/)(4 + 3m + 4i/)(2 + 2i/ + m)(5 + 4i/ + 3m)(3 + 3m + 4i/)a 6 



_ 3 3a 2 9(12i/ 2 + (6m + 12)i/ + 5 + 3m) 2 3(36i/ 2 + 9m 2 + (36/i + 70)i/ + 39m + 40) 

''' ; " '"' _ 16(l + i/) PlP2 + (l + v)(5 + Av + 3^)a 2 Pl 16(1 + v)(2 + m + 2i/)(5 + 3m + 4i/) 



/ 2 20 , 

£3,1 = I- + 6/1 + — t/)a 



3 a 2 3(19 + 21/i + 42i/) 2 3(9m 2 + 12i/ 2 + (24i/ - 9)m - 46i/ - 40) 

Pl + 16 PlP2 + 2(-3 + 3m + 2i/)q 2Pi + (-2 + 3m + 4i/)(-3 + 3m + 2i/) P2 

3(2 + 3m + 6i/)(135m 2 + 21m + 444m^ - 62i/ + 348i/ 2 - 116) 

4(-3 + 3m + 2i/)(m + v)(-2 + 4i/ + 3m)<* 4 Pl 

9(1 + 3m + 6i/)(2 + 3m + 6i/)(135m 2 + 21m + 444mi/ - 62i/ + 348i/ 2 - 116) 

+ 4(1 + 9m + 10i/)(-3 + 3m + 2i/)(m + f)(-2 + 4i/ + 3m)o 4 

/ 22 16 \ 2 

h 6m H v\a , 

V 3 3 / 

3 3a 2 3(12i/ + (61/ + 57)m + 120i/ + 55) 2 3(9m 2 - 12i/ 2 + 12/if - 81m - 210k - 160) 

Pl + 16(9 + i/) PlP2 + 2(9 + f)(M-5)a 2 Pl + 16(9 + y)(/i - 5)(3/i + 2i/ - 14) P2 

9(2 + 3m + 6i/)(48i/ 3 + (96m + 132)i/ 2 + (36m 2 + 792m - 3350)k + 351m 2 - 1347/1 - 1700) 
+ 4(9 + i/)(m - 5) (2i/ + 3m - 6) (3/1 + 2i/ - 14)a 4 Pl 

27(2 + 3m + 6i/)(l + 3m + 6i/)(48i/ 3 + (96/1 + 132)i/ 2 + (36m 2 + 792m - 3350)i/ + 351m 2 - 1347m - 1700) 
+ 4(9 + i/)(m - 5) (9/i + 81/ - ll)(2i/ + 3m - 6)(3/i + 1v - 14)a 4 
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